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ABSTRACT 



From optical spectroscopic measurements we determine that the HD 15407 

Q^ binary system is ~80 Myr old. The primary, HD 15407A (spectral type F5 V), 

exhibits strong mid-infrared excess emission indicative of a recent catastrophic 

^ ■ collision between rocky planetary embryos or planets in its inner planetary sys- 

tern. Synthesis of all known stars with large quantities of dust in their terrestrial 
planet zone indicates that for stars of roughly Solar mass this warm dust phe- 
nomenon occurs at ages between 30 and 100 Myr. In contrast, for stars of a few 
Solar masses, the dominant era of the final assembling of rocky planets occurs 

CN ■ earlier, between 10 and 30 Myr age. The incidence of the warm dust phenomenon, 

when compared against models for the formation of rocky terrestrial-like bodies, 
implies that rocky planet formation in the terrestrial planet zone around Sun-like 

(N) ■ stars is common. 

in 

q ■ Subject headings: circumstellar matter - - infrared: planetary systems — planets 

and satellites: formation — stars: individual (HD 15407A, HD15407B) — stars: 



kinematics and dynamics 



1. Introduction 

Investigating the formation of terrestrial planets outside of the Solar System can give 
insight into the origin of the Earth-Moon system and whether similar planetary system 
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architectures form around other stars. The study of extra-solar terrestrial planet formation 
must currently be carried out through indirect searches for the dusty aftermath of formation 
events. Since the dust is located in the terrestrial planet zone, it will have temperatures of 
~300 K and will emit most of its thermal radiation at ~10 /im. Stars actively undergoing 
terrestrial planet formation will have substantial amounts of warm dust orbiting in their 
terrestrial planet zones and will thus show mid-infrared emission in excess of what one 
would expect from the star alone. Properly establishing the age of such mid-infrared excess 
systems is necessary for comparison with expectations from planet formation theories and 
simulations. 

We were alerted to excess mid-infrared emission around the star HD 15407A through a 
paper entitled "Detection of silica dust in the debris disk around HD 15407" presented at an 
October 2009 Spitzer meeting by H. Fujiwara (E. Becklin, personal communication 2009). 
At this meeting, the age for HD 15407A was suggested to be ~2 Gyr. We felt this age could 
be in error, and initiated our own follow-up observations to independently constrain HD 
15407A's age. Here we report these observations, our revised value for the age of the HD 
15407 binary system, and a discussion of terrestrial planet zone collisions around adolescent 
Solar mass stars. 



2. Observations 

Both components of the HD 15407 system were observed in cloudy condition s at Mauna 



Kea Observatory with the Keck I telescope and HIRES echelle spectrometer (jVogt et al. 



19941 ) . With the HIRES red collimator wavelengths between 4680 and 8700 A were covered 
at a resolution of 50,000 (as determined from the FWHM of single ThAr arclines). The final 
signal-to-noise ratio per pixel, measured at 6700 A, is 350 for HD 15407A and 200 for HD 
15407B. Spectral data were reduced and extracted using both the MAKEE software package 
and standard IRAF tasks and echelle reduction procedures. 



3. Results 

3.1. HD 15407 Stellar Parameters 

Our star of interest, HD 15407A, is part of a visual double star system. Tabled] contains 
information about HD 15407A, the dusty primary star, and HD 15407B, the secondary 
component of the HD 15407 system. For both stars we employ line-depth ratios to roughly 
determine absolute temperature and, for HD 15407B, luminosity. Those interested in the 
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details of determining such parameters are referred to iMelis et all ( 120091 ) 



We deduce that HD 15407A has a temperature similar to that of an F5-type star; its 
broad-band colors (BVIJHK S ) are in agreement with this temperature class determination. 
This suggests that HD 15407A is not significantly reddened by its circumstellar material and 
we adopt a temperature class of F5±l subclass. Spectral analysis suggests that HD 15407B 
is a dwarf star with a temperature class of K3±l subclass and an effective temperature 
of 4600 K. Broad-band colors for HD 15407B instead suggest that this dwarf star has an 
earlier spectral type of K0-K1. We adopt a spectral type of K2 V with an uncertainty of ±1 
temperature subclass. 



3.2. Excess Infrared Emission from HD 15407A 

Excess infrared emission is d etected towards HD 15407A by IRAS in the 12 and 25 
(Urn bands (jOudmaijer et q/.lll992l ). Archival Spitzer IRS spectroscopy of HD 15407A (AOR 
26122496, PI Fujiwara) confirms the excess (see Figured]). We note the presence of solid 
state dust emission features in the IRS spectrum. Figure [2] compares the IRS spectrum of 
HD 15407A with those of other stars that host substantial quantities of terrestrial planet 
zone dust. We defer quantitative analysis of these spectral features to later publications. 

To estimate the dust temperature and the fraction of the stellar luminosity reradiated 
by the dust (r=LiR./Lboi) we fit optical an d near-infrared measur ements out to K s -band with 
a synthetic stellar atmosphere spectrum (IHauschildt et q/.lll999l ) and a blackbody to model 
the dust excess (see Figured]). In this way we find that the fraction of stellar light reradi- 
ated by the circumstellar dust is ~5.7xl0 -3 and that the grains radiate with a blackbody 
temperature of ~500 K. Such grains will reside at a distance of ~0.8 AU from HD 15407A. 

We note that both the IRS spectrum and an IRAS upper limit at 60 /im of 347 mJy 
(4cr) suggest that there is no cold dust orbiting HD 15407A. 



3.3. Age of the HD 15407 System 

To determine the age of the binary system we obtained and analyzed high-resolution 
echelle spectra for both stars and searched for X-ray detections in the literature. From these 
data we estimate the system age from the lithium content in the stellar photospheres, velocity 
widths of a bsorption lines, Galactic sp ace motion, and chromospheric activity; details can 
be found in IZuckerman fc Song! (120041 ) . The lithium 6710 A absorption feature is strongly 
detected in the spectra of both stars (Tableland Figure [3]). Lithium content in a stellar 
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atmosphere is mainly determined by the star's age and mass. These strong detections are 
consistent with roughly Solar mass stars of ~30-100 Myr age (see Figure[3]). We note that the 
lithium equivalent width (EW) measured for HD 15407A (~94 mA) is significantly stronger 
than the range of Li I EWs (~40-70 mA) quoted for similar spectral type Pleiads in Figure 
3; this suggests that the HD 15407 system may be younger than Pleiades cluster stars. 

The ROSAT All-Sky Survey (RASS) detected X-ray emission from HD 1 5407B. Analysis 
of thi s X-ray emission (TableEDand F i gureBJ) suggests a star of ~100 Myr age (IZuckerman fc Song 



2004ft . We note that ISuchkov et all (1200a ) claim a detection of HD 15407A by the RASS. 



Upon further investigation, it is revealed that the putative detection has an angular separa- 
tion of >5' from HD 15407A, well outside the RASS positional uncertainties. We conclude 
that HD 15407A was not detected by ROSAT. We estimate an upper limit for its fractional 
X-ray luminosity in the following manner: HD 15407B is detected by ROSAT with ~4a 
significance. We take this level of X-ray flux as our upper limit for HD 15407A. Thus, we 
derive the upper limit to HD 15407A's fractional X-ray luminosity by scaling HD 15407B's 
fractional X-ray luminosity by the ratio of the two star's total luminosities. This upper limit 
is reported in Table [T] and plotted in Figure H] where it is shown to be consistent with stars 
of ~100 Myr age. 

Velocities of the HD 15407 system toward the center of our Milky Way galaxy, around 
the Galactic Center, and perpendicular to the Galactic plane (U, V, W) are ca lculated from 
Table [T] sky positions and proper motions, the Hipparcos measured distance (Ivan Leeuwen 



20071 . where we assume HD 15407A and B are bound as is suggested by their common proper 



motions and the optical echelle measured radial velocities, see Table [I]), and heliocentric 
radial velocities. UVW is computed for each component of the HD 15407 system and these 
two results are averaged together to produce the values reported in Table [TJ Comparison of 
these com puted UVW space motions to those of known nearby stellar associations (see Figure 
E]and e.g. . IZuckerman &: Song||2004l ; iTorres et a/.ll2008l . and references therein) suggests that 
the HD 15407 system is a member of the AB Dor moving group (UVW~— 8, —27, —14 
km s _1 ). The distance to the HD 15407 system and its position in the plane of the sky 
are consistent with it be ing a member of the AB Dor moving group (e.g., see Table 5 of 
Zuckerman fc Songj|2004ft. Age estimates for the AB Dor moving gr oup range from ~50-120 



Myr (IZuckerman et al. 



2001 



Luhman et aZ.ll2005l ; ITorres et a/.ll2008l . and references therein). 



Based on the above analysis, we suggest a best estimate age for the HD 15407 system 
of ~80 Myr with a conservative range of plausible ages from 60-120 Myr. We note that the 
measured vsmi for HD 15407Elj, ~4 km s _1 , is surprisingly low for a star of such age and 



1 vsini for both stars was measured from the full-width at half-maximum depth (FHWM) of single ab- 
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spectral type (IZuckerman Sz Song||2004l . and references therein); it is possible that we view 
HD 15407B in a pole-on geometry. 



4. Discussion 

Here we use our age estimate for the HD 15407 system to help constrain the era of rocky 
planet formation in the terrestrial planet zone around stars of Solar mass. We consider stars 
with fractional infrared luminosities r>0.01 as very dusty and likely to have undergone a 
recent collision between two planet a ry em bryos or planets such as that postulated to have 
formed Earth's Moon. iRhee et al\ (120081 ) considered this problem, but included the very 
dusty G-type star BD+20 307 in their analysis. Subsequent to their paper, it was discovered 
that BD+20 307 is an old main seq uence star and hence its circumstellar du s t is not re 



l ated to terrestrial planet formation (jWeinbergerl l2008t IZuckerman et q/.ll2008l ). IRhee et al. 



( 120071 ) and C. Melis et al. (2010, in preparation) have investigated "intermediate mass" stars 
(spectral types late-B through early-F and masses from ~8 to ~1.5 M Q ) with large masses 
of warm dust in their terrestrial planet zones. The half dozen or so very dusty intermediate 
mass stars currently known have ages between 10 and 30 Myr. This suggests that for such 
stars the era of most intense accretion of material onto rocky planets is over by about 30 
Myr. Such a conclusion is similar to that drawn by ICurrie et all (j2008f ) who considered an 
ensemble of young intermediate mass stars but with only modest amounts of warm dust. We 
now show that the corresponding era for solar-mass stars, by contrast, lasts for ~100 Myr. 

Based on IRAS and Spitzer surveys we are aware of four adolescent Solar-typ e stars with 
t>0. 01: AB Dor member H P 15407A (this pap er); Pleiades member HD 235 14 (IRhee et al. 
2008 ): M47 member PI 121 (JGorlova et a/.ll2004h : and NGC 2547 member ID8 dGorlova et al. 
20071 ). The ages of these stars range between 35 Myr (NGC 2547) and about 100 Myr (the 
other 3 stars). Interestingly, two of the four young, very dusty Sun-like stars mentioned above 
are in wide binary systems (HD 23514, D. Rodriguez et al. 2010 in preparation; HD 15407, 
this work). The other two systems (ID8 and PI 121) have not been exhaustively searched 
for binary companions; confirmation of any such binarity would be suggestive (although 
still in the realm of statistics of small numbers) of a link between binarity and catastrophic 
terrestrial planet collisions. 

We note the relative absence of similar numbers of warm dust stars with ages from 



sorption lines in the Keck HIRES spectra; the intrinsic FWHM resolution for the HIRES setup used was 
^6 km s , a value we subtract in quadrature f rom the FWHM measured in the spectra when determining 
vsini (see Equation 6 in IStrassmeier et ali\199(jt ). 



-6- 



200 to 800 Myr. No stars i n four stellar associa tions with age s in this range — the H yades 



fjCieza et al 



20081). Praesepe (IGaspar et a/.ll2009l ). Carina- Near (jZuckerman et al\\200w , and 



Ursa Major (jKing et a/.ll2003l ) — have substantial quantities of warm dust in their terrestrial 
planet zone. Together these four groups contain ~200 Solar-type stars. In the field, IRAS 
has surveyed many more stars with ages of ~200 to 800 Myr compared to the 30 to 100 Myr 
range, yet none in the former age range have been identified with luminous warm dust. 

We estimate the frequency of catastrophic accretion events during rocky planet forma- 
tion in the terrestrial planet zone by considering the detection rate of very dusty stars in 
the field and in young clusters. Very dusty nearby field stars were better investigated in th e 
all-sky IRAS survey than with Spitzer's pointed capabilities. As noted by lRhee et all ( 120081 ). 
IRAS could have detected a main sequence dwarf with r>0.01 and spectral type between 
F4 and K0 out to about 150 pc. There are 400 or so such field stars with ages ~100 Myr 
that IRAS could have detected, but of these, only HD 15407A has sufficient dust to actually 
have been detected by IRAS. Simila rly, among the youn g clusters st udied with Spitzer 
NGC 2232J 



Currie et all 



2008)), M47 (IGorlova et q/.ll2004h. NGC 2547 (IGorlova et al. 



2003) 



IC 2391 (ISiegler et o/.l 120071 ). and the Pleiades flStauffer et q/.ll2005l : IGorlova et al\\200& ] 
the three very dusty cluster stars listed previously represent about one Solar- type star in 
200 (roughly consistent with the determination for Solar-type cluster stars by iBalog et al 



20091) 



Thus, the high r warm dust phenomenon manifests itself at about one adolescent Solar- 
type star in 300. If all F4-K0 stars display this phenomenon as adolescents, then the lifetime 
of the phenomenon at a typical Sun-like star is about 300,000 years. 

To understand this frequency of th e warm dust pheno menon we consi der a dust produc- 
tion model similar to that described bv iRhee et all (boost ) for HD 23514. iRhee et all (boosh 
analyzed existing observations through October 2007 with the aid of a model of colliding 
planetary embryos due to C. Agnor and E. Asphaug (jAgnor et ctZ.I Il999l : lAgnor fc Asphaug 
20041 ; lAsphaug et al\\200w and earlier researchers. The idea is that following a catastrophic 
collision of two rocky planetary embryos or planets, fra g mente d debris covering a range of 
sizes goes into orbit around a central star. IRhee et al\ (120081 ) assumed that collisions be- 
tween debris would be frequent enough to establish an equilibrium size distribution where 
the number of particles of radius "a" is proportional to a -3 ' 5 . In this situation, both the 
mass carried by particles of radius a and the collision time are proportional to the square 
root of a. Most of the mass is carried by large fragments while small particles are responsible 
for most of the excess infrared luminosity and collide the fastest. The radius of the smallest 
dust particles in orbit around a very dusty star is set by radiative blowout and, for HD 
15407A, will be ~1 /im. The collision time for such particles will be about 1 year divided 
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by 10xLiR/Lboi or ~10 years. If the largest initial fragments of a catastrophic collision have 
radii ~100 m, then their collisional lifetime is ~100,000 years which is within a factor of 3 of 
the lifetime (estimated in the previous paragraph) required if all Solar-type stars are form- 
ing rocky planets in their terrestrial planet zone. Should the formation of rocky terrestrial 
planets involve more than a single catastrophic collision, then the model-estimated warm 
du st phenomenon life time could be more than 100,000 years at each Sun-like star. As noted 



by iRhee et al\ ( 120081 ). the total mass lost over this interval of time — due to a collisional 



cascade followed by radiative blowout — will be about the mass of Earth's Moon. 



5. Conclusions 

Based on Keck HIRES optical spectroscopy, X-ray flux, and Galactic UVW space mo- 
tions we estimate the age of the HD 15407 binary system to be ~80 Myr and identify the 
system as likely belonging to the AB Dor moving group. Collecting the full sample of known 
very dusty debris-disk hosting stars, we find that the era of final catastrophic mass accretion 
for rocky planets around stars of roughly Solar mass occurs later than around intermediate 
mass stars. Provided that rocky planet formation in the terrestrial planet zone of most So- 
lar mass stars is accompanied by at least one and possibly a few catastrophic collisions of 
planet-size objects within a 100 Myr interval, then the detection rate and age range of very 
dusty Solar-type stars — observed to be about one in 300 for ages between 30 and 100 Myr 
— can be accounted for. Our conclusion is that existing observations of such dusty stars 
indicates that most Sun-like stars are accompanied by rocky planets in the terrestrial planet 
zone and that collisions among such forming planets will occur, as is expected in theoretical 
simulations of planetary origins. 

We are indebted to George Herbig for obtaining the HIRES data for us. We thank 
Eric Becklin for pointing out this stellar system to us. Some of the data presented herein 
were obtained at the W.M. Keck Observatory, which is operated as a scientific partnership 
among the California Institute of Technology, the University of California and the National 
Aeronautics and Space Administration. The Observatory was made possible by the generous 
financial support of the W.M. Keck Foundation. This research has made use of the VizieR 
service and of data products from the Two Micron All Sky Survey. This research was 
supported in part by a NASA grant to UCLA. 
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Fig. 1. — Significant mid-infrared excess emis sion is detected towards HD 15407A. The solid 
brown curve is a synthetic stellar spectrum (IHauschildt et q/.lll999f ) for a 6500 K effective 
temperature star. The brown dotted curve represents a dust continuum emission blackbody 
fit to a dust temperature of 500 K. The black dashed line, the sum of the above two curves, 
is a reasonable fit to all data points. The total luminosity of the excess, determined by 
integrating under the data points between 1 /im and 100 /zm, is ~0.6% of the luminosity of 
the star. The BV data points are from the Tycho-2 catalog. The green data points are from 
the IRAS Faint Source Reject Catalog and have not been color corrected. The horizontal 
bars indicate the filter bandwidths. Overplotted is the reduced, extracted, calibrated Spitzer 
IRS spectrum. 
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Fig. 2. — Top Panel: Comparison of HD 15407A's IRS spectrum to those of othe r stars that 
host 9 /iin spectral features. HD 15407A (r~0.6%) and HD 23514A (r~2%; see JKhee et al. 
20081 . D. Rodriguez et al. 2010 submitted) are ~100 Myr old mid-F type stars while HD 
172555 (r~0.08%) is a ~12 Myr old late-A type star. HD 172555's IRS spectrum has been 
analyzed in great detail and is dominated by features from amorph ous olivine and pyroxene, 
crystalline pyroxene, and silica ( iChen et q/.ll2006l ; iLisse et aZ.ll2009h : the latter dust species is 
suggested to result from high velocity impacts. Bottom Panel: I RS spectra of other h ot-dust 
hosting stars. ID8 (r~2%) is a ~35 Myr old e arlv-G type star (JGorlova et d.l l200~7l). P1121 
(r~l%) is a ~100 Myr old late-F type star (JGorlova et a i]|2004h, BD+ 20 307 fr~4%) is 
a >1 Gyr old spectroscopic binary composed of two early-G type stars (jZuckerman et al. 



20081 ). HD 113766A (r~1.5%) is a 

EF Cha (r~0.1%) is a ~10 Myr old late-A type star f lRhee et al. 



16 Myr old early-F type star (IChen et aZ.1 120061 ) . and 

20071 ) . HD 113766A's 



IRS spectrum has been analyzed in great detail a nd is dominated by amorphous olivi ne and 
pyroxene features and crystalline olivine features (IChen et q/.ll2006t ILisse et aZ.ll2008l ) . 
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Fig. 3. — Top Left Panel: Keck HIRES spectra obtained on UT 29 Nov 2009 and covering the 
Li I A6710 (vacuum) region for HD 15407A and B. Strong lithium lines are detected in both 
stars (see Tabled]). These spectra are continuum normalized and offset by arbitrary values. 
Wavelengths in this figure are plotted in the heliocentric reference frame and are in vacuum. 
Bottom Left Panel: Comparison of the measured Li I EWs of HD 1540 7A and B (Table [1]) 
to th ose of moving groups with similar age (adapted from Figure 3 of IZuckerman Sz Song 
20041 ). The Pleiades age is ~100 Myr, the Tuca na association age is ~3 Myr, and the 
AB Dor association age is ~80 Myr (see text and IZuckerman fc Songll2004l ). Right Panels: 
Comparison of the Galactic UVW space motion of the HD 15407 system (red star symbol) 
with the UVW space motions of (1) the cluster o f plus symbols nea r the c enter of the figures 
that are AB Dor moving group members from IZuckerman et all (120041 ) and (2) triangle s 



that are white dwarfs with mean total ages comparable to that of the Sun (IZuckerman et al. 



as defined in 



20031 ; iBergeron et a/.ll200lh. The rectan gles represent the "good" UVW box for young stars 



Zuckerman fc Song! (120041 ) . 
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Fig. 4. — Comparison of the calculated fractional X-ray luminosities of HP 15407A and B 
(Tab le [1]) to those of two clusters of known age (adapted from Figure 4 of IZuckerman fc Song 
20041 ). HD 15407A was not detected by the ROSAT satellite, the plotted value is an estimated 
upper limit (see text). The Pleiades age is ~100 Myr while the Hyades age is ~600 Myr. 
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Table 1: Parameters of the HD 15407 binary star system 



Individual Stars 






HD 15407A 


HD 15407B 


RA (J2000) 




02 30 50.65 


02 30 48.58 


DEC (J2000) 




+55 32 54.4 


+55 33 06.4 


Sp. Type 




F5 V 


K2 V 


Vmag 




6.95 


9.58 


T cff (K) 




6500±100 


4830+150 


pmRA (mas yr~ 


V ) 


+80.6 ± 1.0 


+83.1 ± 1.0 


pmDE (mas yr~ 


V ) 


-96.7 + 1.1 


-93.6 ± 1.1 


RV (km s" 1 ) 




-10.44 ± 0.39 


-10.02 ± 0.17 


Li I A6710 EW (mA) 


94+4 


174+5 


vsim (km s" 1 ) 




20 


4 


Lx/Lbol 




<1Q -5.0 


10 - 4 .i 


System 



Distance from Earth (pc) 

Separation Between Stars 

UVW Space Motions (km s" 1 ) 

Age (Myr) 



54.7+2.1 

21.25" (projected separation of 1160 AU) 

-11.5, -28.8, -14.0 



so; 



-10 

-20 



Note. -- Epoch J2000 positions are from the 2MASS catalog ([Cutri et al\\2003i) . The highest fidelity 



proper motions come from th e Tycho-2 catalog ( Hog et gZ.I l2000h for HD 15407A and from the UCAC3 



catalog ( Zacharias et aZ.ll2010f ) for HD 15407B. Radial velocities (RV) for both stars were measured by cross- 
correlating our HIRES spectra with HIRES spectra of a star from Table 2 of lNidever et all (|2002h . See the 
text for discussion of other stellar parameters. 



